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1. Introduction 


On atom chips, coherent atomic waveguides are important components needed for the construction 
of certain types of matter wave interferometers that use cold magnetic atoms. Results concerning 
the use of the time reversal symmetry (Tinkham, 1964) to solve the quantum behavior of a spin one- 
half particle guided in a type of atom waveguide often used on atom chips are presented here. 

The Hamiltonian for the atom guide is symmetric under time inversion, if the reversal of the 
external magnetic fields is explicitly included in the reversal process. As described in a previous 
report (Golding, 2010a), the Frobenius power series method can be used to derive a pair of power 
series solutions that are well behaved at the coordinate origin. These two solutions are related to one 
another by the time reversal symmetry of the system. Two other solutions that are not finite at the 
origin are rejected for physical reasons in this waveguide problem. 

These series solutions are useful as they allow the basic symmetry properties of the transverse 
modes of the guided atoms to be studied analytically. In addition, these series can be readily 
summed for a large range of radial distances using arbitrary precision arithmetic techniques so that 
the behavior of guided modes can be studied numerically. 

In order to form eigenstates and calculate eigenvalues, the pair of physically acceptable solutions 
derived in this way must be combined to satisfy some type of boundary condition relevant to the 
problem of interest. Others have essentially used an open boundary condition to describe a 
quadrupole guide of infinite transverse extent. This leads to an ambiguity in the eigenstates that may 
be important in understanding the behavior of these guides in practice. In our current work, a hard- 
wall boundary condition is used, establishing a closed boundary. This forces both components of 
the spinor solutions to be zero at some fixed radius representing the boundary. When hard-wall 
boundary conditions are used, a complete system of orthogonal eigenstates is found that satisfies the 
Sturm-Liouville conditions. Further work is needed to understand the relationship between these 
two sets of boundary conditions as the limit of a large radius for the hard wall is taken. 

Since explicit series solutions are developed, an ordinary differential equation (ODE) solver is not 
needed to calculate the wavefunctions or eigenvalues. This is extremely useful as the fourth order 
differential equations involved are stiff and difficult to handle with nonnal ODE techniques. Both 
the energy spectra and the wavefunctions can be calculated directly using the series technique for 
various regimes of magnetic field strength, allowing the direct study of the assumed adiabatic 
behavior (Sukumar,1997) as the longitudinal bias field (Ioffe field) is increased from zero to 
relatively large values. 
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2. The Atom Guide Model 


To study the dynamics of an atom in a magnetic waveguide, the Schrodinger equation describing a 
single neutral magnetic atom trapped in the magnetic guiding field must be solved. To minimize the 
initial complexities, the initial model is chosen to be a neutral spin Vi atomic level. Two alkali 
atoms, 40 K and 6 Li, are fermions and have a spin !4 level. The important class of bosonic alkali 
atoms will be treated by a similar technique in a later work. 

3. Hamiltonian for a Magnetic Atom in a Magnetic Waveguide 

The quantum description of an atom moving in a spatially dependent magnetic field must include 
both internal and external degrees of freedom. This means that the momentum, position, and spin 
degrees of freedom must all be treated as quantum operators. Using m for the mass of the atom and 

p 1 

M for the magnetic moment, the Hamiltonian is written as the sum of the kinetic energy, ~ and 

the interaction energy of the magnetic moment and the field B . The result is 

2 

H = —— M B(r) 

2 m \ (1) 

The magnetic field B(f ) is independent of the z-direction and is taken to be an ideal quadrupole 

field equation (equation 2) extending to infinity in the x-y plane and uniform along z. An additional 
uniform bias field B 0 is added in the z-direction to help control possible spin-dependent losses that 
may occur at the zero field point at the center of the guide. The spatial dependence of the ideal 
quadrupole field is given by 

B(r) = B l (-xx + yy) + B 0 z ^ 

In the simplified model of the atom described by equation 1, the magnetic moment of the atom is 
just the magnetic moment of a single outer electron. The atomic model used here can be thought of 
as an alkali atom with a spin zero nucleus and a total mass m. The magnetic moment is M = yS , 
where y is the gyromagnetic ratio of the level considered and S is the spin angular momentum of 
the atom. Since a spin Vi system is being considered, the spin angular momentum is proportional to 
- fi 

the Pauli matrices, S = — a . 

2 

The Schrodinger equation for the guided atom eigenstate ^ E is 
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( 3 ) 


Zm 

Equation 3 is made dimensionless by choosing a length scale such that x —> Xx, where the new x is 

n 2 k 2 

2m 

where k = — is the wavenumber of the optical field used to cool the atom. This is a rather arbitrary 

X 

choice at this stage, but the recoil energy is a convenient scale for atoms that have been laser 
cooled. This is done after dividing through by the recoil energy, pulling out the relevant factors, and 

_ kB x X 

defining the dimensionless transverse field parameter U ~ 7 ~ , the longitudinal field parameter 

lb R 

_ HB 0 _ E 

b o ~ 7 -r~r , and the scaled energy £ - — . 

Zh R E r 


4. Time Reversal Symmetry of the Hamiltonian 


Pr 2 - _ frk 2 

dimensionless. Then — ^ ^ e ( r ) “^ TZ ^ ^ ( r ) and the recoil energy is defined as E R 


2m 


2m 


The Hamiltonian above (equation 1) is invariant to the operation of time reversal if the external 
field B is reversed as well. This is simply because when the spin, S , reverses sign under time 
reversal, B must reverse sign to leave the Hamiltonian unchanged. The usual time reversal 
techniques do not act on an external magnetic field. An operator M is defined to explicitly reverse 
the external magnetic fields in the calculations, thus MBM 1 = —B . When the dimensionless field 
parameters are present in an expression, the field reversal operator simply changes the signs of both 
b 0 and b x . This operator is then combined with the standard time reversal operator, T = e ' S ''' T ' K , 
where K is the conjugate operator used in scalar time reversal problems, to form a symmetry 
operation of the system (equation 1). The complete time reversal symmetry of the Hamiltonian is 
then expressed as 

(MT)H(MT ) 1 =H . (4) 

Using this operator and any given initial state u) an additional state I ¥ ’) called the time reversed 
state can be fonned as | ¥ ’) = MT \ y) _ 
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5. Eigenstates Common to H and A z 


As described in Golding (2010a), the equations for the radial wave functions of the quadrupole 
guide are derived by making use of the angular symmetry of the quadrupole magnetic field. This 
symmetry is expressed as the conservation of the angular momentum, A _=L_- S, (Hinds, 2000 

and 2001; Lesanovsky, 2004; Bergeman, 1989; Golding, 2009). Our approach is to find 
eigenfunctions that are common to both the Hamiltonian, H and the angular momentum, A _. This 

approach effectively separates the angular and radial dependence of the problem, leaving only a 
coupled system of radial equations (equation 6) to solve. These equations are solved exactly using 
both series techniques and the time reversal symmetry of the Hamiltonian (equation 1) to develop 
the general solution of the Schrodinger equation. The general solution must then be made to satisfy 
appropriate boundary conditions so that eigenvalues and eigenfunctions can be calculated. 


The component fonn defining the radial spinor components for a state with angular momentum pi is 




R](p)e ,,p/2 
R u {p) e 1,1,11 


When the assumed fonn (equation 5) is used in the Schrodinger equation (equation 3) the angular 
dependence drops out leaving only the following radial equations to be solved: 


d;R + + -d p R + - 

P 


(p + l/2f 


R + +(z + b 0 )R + = b x p R 


d;R+-d p R- 

P 


(A-1/2 f 


R + (s - b 0 )R_ = b x pR + . 


In the case of a spin half particle with integral orbital angular momentum, the pair of equations must 

l 3 5 

be solved for the allowed half-integral values of the angular momentum ^ ~ ± 2 ,± ^. ,± 2"' ■ 


6. Frobenius Series Approach 


When the system of radial equations (equation 6) is uncoupled, a fourth order differential equation 
results that has a regular singularity at the origin. This equations can be solved using the Frobenius 
power series technique (Ince, 1956; Bender, 1978). The equation for R + is 
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4 d 4 _ 

P ~r^ R + + 


dp 

( s A 

2 p'e- — + 2p 2 

\2 J 


P 


d 2 


dp 


R + 


2 + 


d 


-2p 2 b Q +(5 + 4p 2 +6p}pj-j— R + + 

~b\ p b + {s 2 -b 2 )p A +ib Q - 2 p 2 s + 2pb 0 + ^s p 2 -y // 2 -9//-^ + // 4 


i? =0 


( 7 ) 


A similar equation for R can be found simply by replacing R with R and changing the signs of 
both p and bo. This is readily understood by inspecting the behavior of the coupled radial equations 
(equation 6) under this transformation and is due to the time reversal symmetry discussed 
previously. 


The Frobenius series technique is applied to equation 7 by assuming a modified power series 
solution for R+ of the form 


*-° , (8) 

where the coefficients c„ and the index a are u nkn own except that co = 1. This is substituted into 
equation 7 and the resulting system is solved resulting in the following recurrence relations for c„, 


c n-b b x + c n -4 (V -)- C n _ 2 (p + 5/2-n-a)(2sp-2b 0 -5s + 2sn + 2sa) 

c —-----. (9) 

(2// + 1 + 2/7 + 2<t)( 2//-7 + 2/7+ 2cr)(2// + 5 -2«-2cr)(2// + l-2«-2cr) 

For n = 0, the denominator of equation 9 is the indicial equation that is used to determine the 
allowed values of a in equation 8. The necessary conditions are just that co = 1 and c n = 0 for 
n < 0. The corresponding component R is found by substituting equation 8 into the first line of 
equation 6 and solving for /?_ producing 

R =I >.{(«+Ad/' +1/2 )V +£ f^ 1 '+ < 10 ) 

n =o ^ ' b x p 

Since the allowed values of a differ by integers, the maximal index is the only index that is 
guaranteed to produce a solution that is finite at the origin (Ince, 1956). This solution is called a 
type one solution and the solution obtained from this by time reversal is called a type two solution. 
Type one and type two functions are denoted with a leading superscript such as V or R in the 
following. 
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In order to use equation 9 over the full range of angular momenta, pi, the fact that the maximal index 
cr is a function of pi must be considered. The maximal index a is defined as a function of pi by 


cr 




pi> 0 
pi < 0 


(ID 


7. Time Reversal and the Two Solutions 


Equations 8 through 11 define a single solution to equation 6 for the specific angular momentum pi 
that is called the type 1 solution, V /( . The type 2 solution, V"“, is derived directly from the type 1 
solution using the time reversal symmetry of the system. In terms of the radial wavefunctions, the 
two solutions are defined as 


r 

r 


PR»(b t ,b,)e»'\ 
' 2 K(b < ,,b,)e‘ ,n ' 

pR"(b t ,b t y rn / 


and 


( 12 ) 


The full time reversal symmetry MT acts on )//■“ to give V“, but the fields and angular momenta 
have to be properly reversed as in 

( Inul L L\.>/2\ f 2n -11 (u 


MT\^ =e i ^ 


where 


~'Ra(-b,-by 

k 'r:(-k ,-by vl , 




y (byy 


2 R7 (b„,b,)e- 


'<p! 2 


= v- 


(13) 


MT = -icjMK = M 


0 -1 

1 0 


K. 


Thus, the time reversal operator creates a second spinor solution independent of the first as can be 
demonstrated by calculating the Wronskian of the solution pair in equation 13. The action of the 
time reversal operator on the spinors of different type is summarized as 


MT V" (/>„./>,) -T >")- V" (b,A ) and 
MT V (b 0 ,b,)=T V" = - V (KA)- 


By inspection of equation 13 the component radial functions of the two solution types are related by 

'R:(b,Ay 1 R: u (-b„-b l ) 

'R'!(b„A) = - 2 RA(-b a ,-b t )' 
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Application of equation 15 to the type one series solutions gives the series solution for the other 
radial components in terms of the single set of coefficients c „ of equation 9. These series are 
displayed in equation 16: 


X K (b 0 ,b x ) = p <M) Y,c„ (■ cr ( n ), mA ,b x ) p n 
P 


n =0 
o-(//)-3 oo 


' R -( b oA) = iL T —(cr(//),//,h 0 ,h 1 )|(« + a(//))" -[p + H lf\p n + —p~p- X R+ {b 0 , b i) 


6, 


1 n =0 


V 


( 16 ) 


oo 

2R - { b oA) = Yj C n( CJ (-p)-P ~ b 0 ’ ~ b i ) P" 


2 K(boA)- 


p 


n=0 

oo 


b 


1 n=0 


'E c A a (-p)’-P’- b 0’- b i ){(« + -(//-l/2)“ (/?" 


s-b. 


b i P 




The series expansions are used in the general form for the spinors (equation 12) to produce a pair of 
basic solutions to the Schrodinger equation (equation 3). 


8. Eigenfunctions and Energies 


The coefficients c n (equation 9) are functions of the energy s . Some type of boundary condition 
must be used with the radial functions (equation 16) to detennine allowed values of the energy. One 
possibility is to use a hard-wall boundary condition and force the general spinor to be equal to zero 
at some radius. Another possibility is to assume there is no hard boundary and to set the general 
solution equal to the asymptotic solution at some convenient distance from the guide center. 
Experimental implementations of a quadrupole guide will not extend to infinity and some 
understanding or model of the physical boundary condition will be needed in order to adequately 
model the actual guide behavior. 

The simplest type of boundary condition is the hard-wall condition. This condition requires setting 
the general solution equal to zero at the wall. To satisfy this boundary condition, a linear 
combination of type one and type two solutions is fonned and eigenenergies are found such that 
both the general form for both the plus and minus components of the spinor are simultaneously 
zero. 

The general linear combination of the two types of solution can be fonned using a single parameter 
9 since there will eventually be a normalization requirement on the states. Using this parameter, the 
general solution is 
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( 17 ) 


y/ M (i 9 ) = cos(<9/2) V" + sin (0/2) V 

f (cos(6'/2) 1 <+sin(^/2) 2 <)e' W2 " 

[(cos(0/2) +sin(6>/2) 2 /A)a^ 

Assuming the boundary condition is independent of V , there are two equations that must be solved 
simultaneously for a common value of 0 and £ . Using A, as the radius of the hard wall, the 
equations detennining the eigenvalues are 

cos(6>/2)X' , (ff,A)+sin(6>/2) ! ^'(s, A ) = 0 

(1 o) 

cos(0/2) 1 R M (£,p o ) + sm(0/ 2) 2 R fl (£,p 0 ) = 0 

When these equations (equation 18) are solved a set of orthogonal spinor radial wavefunctions are 
found that may be used for modeling the behavior of atom guides under hard-wall boundary 
conditions. To model open boundary conditions of infinite extent, a similar procedure must be 
performed using the asymptotic solutions of the system (Golding, 2010b). 

When the Ioffe field is large and the boundary conditions are open, discrete energy levels may exist 
but the transverse modes are not well defined. The lack of control or consistency of the transverse 
mode may lead to difficulties in the construction and operation of beam splitters and other matter 
wave components needed for guided matter wave interferometry and may be one reason that single 
mode operation in these guides has not yet been convincingly demonstrated. One solution to this 
problem would be to develop a new waveguide design that includes both a magnetic quadrupole 
field and some type of conservative hard-wall boundary defined at fixed distance from the guide 
center. This may be possible using a combination of optical fields and miniature magnetic 
quadrupole fields, for example. 


9. Conclusions 


The multiple solutions found in other atom waveguide work can be easily derived using the time 
reversal symmetry of a waveguide system. By organizing the series solutions using time reversal, a 
more complete understanding of guide behavior can be obtained. For example, at large Ioffe fields, 
the two solutions correspond physically to a bound mode and an unbound mode. The bound mode is 
a quantum state that is essentially attracted to the center and its time reversed partner is an unbound 
mode that is repelled from the center. These two modes must be combined to meet prescribed 
boundary conditions resulting in a specific mixture of bound and unbound modes that can be 
directly related to the lifetime of waveguide states. 

This symmetry will be useful in more complicated devices currently being developed, such as 
magnetic beam splitters. For example, if one leg of a beam splitter favors the propagation of a 


8 





bound state and the other leg favors an unbound state, the splitter will not work as desired. It is also 
hoped that the time reversal technique presented here will be useful in the analysis of atom guides 
employing higher spin systems. In such systems, there will be many possible modes involving 
different spin states. Each state will have a time reversed partner that can easily be obtained, 
reducing the need to develop many different independent series solutions. Time reversal techniques 
reduce the analytic and numeric work needed to obtain a general solution as well as enhances the 
conceptual understanding of atom guide systems. 
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